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We investigate the resonance steps, spatiotemporal dynamics, and dy- 
namical phase diagrams of the dc-driven diatomic Frenkel-Kontorova model. 
The complete resonance velocity spectrum is given. The diatomic effects re- 
sult in each resonant state being characterized by two integer pairs, {ki,k2) 
and {ki,k'2)- In the high- velocity region the linear response of the average 
velocity , v, of the chain to the driving force F is often punctuated by subhar- 
monic resonances {ki > k2). There are two kinds of nonlinear response regions 
in the high-velocity region. A new physical interpretation of the mean-field 

treatment is presented. 
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The static and the transport properties of Prenkel-Kontorova(FK) [1] systems have been 
subjects of intensive studies [2-4]. In these studies, one of the most important parameters 
is the winding number uj [5] . Starting from uniform shding states, Strunz and Elmer [6] pre- 
sented the main, superharmonic resonant conditions while Zheng et al [7] gave the complete 
resonance condition. The main studies in Refs. [6] and [7] are limited to the region of low 
velocity and to the case where the chain is composed of the same atoms. 

In this paper, we study the dc-driven FK system where the chain is composed alterna- 
tively of two different kinds of atoms. The response of the average velocity of the system 
to the driving force, the resonant steps, and the spatiotemporal dynamics are studied ana- 
lytically and numerically. Our studies shed some different light on the dynamical behavior 
of the driven FK systems because (i) diatomic effects which result in new resonance modes 
are shown; (ii) regions far away from those addressed in previous investigations are treated; 
and (iii) a new physical interpretation of the mean-field treatment is presented. 

We consider a uniform, quasi-static force field F{t) acting on a FK system in the presence 
of a damping proportional to the velocity and where the atoms can exchange energy with a 
thermal bath. When the FK chain is composed of more than one kind of atoms, the motion 
equation of the system may be expressed as 

rrij Xj — —rrijr] Xj +{xj-i — 2xj + Xj+i) 

-^sm{27rXj)+F . (1) 

In this paper, we choose mi = nis = ••• = !, 1712 = "024 = ■ ■ ■ = m where m > 1. The ground 
states of this diatomic FK model have been studied in one of our previous papers [8] .The 
diatomic effect results in a second critical point. 

When K = 0, there is no nonlinear interaction in the system, and the system mobility 
shows a linear response B — 1/r]. When K > 0, the interactions between the atoms and 
the external potential are nonlinear. The commensurate ground states are always pinned. 
For an incommensurate ground state with the winding number uj, there exists a critical 
value Kc{uj) at which the Aubry transition occurs [9]. When K < Kc{u!), the system is 
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unpinned; when K > Kcicu), the system is pinned. For the latter case, when F is less 
than the depinning force F^, the system is in the locked state. When F > K/{2t^), the 
system should behave approximately as a homogeneous one, and B{F) reaches the final 
value — l/rj (linear response). For an intermediate force, the mobility depends on the 
driving strength (nonlinear response). 

The average sliding velocity of the chain reads 



For an incommensurate structure with winding number cu, when K < Kc{oj) at which the 



dissipation, the chain may slide without any applied force. If the particles are charged, this 
state corresponds to the superconductivity state. 

For a numerical solution to Eq. (1), we use the standard fourth-order Runge-Kutta 
method. Periodic boundary conditions, xat+i = x\-\- M, xq = xn — M, Xo=X]sf, X]sf_^i=Xi, 
are imposed, where M, N are positive integers and N is an even number. In the simulation, 
we start with the ground state and F = 0. We increase the dc force F by small steps AF 
up to a value Fmax- Then, we decrease it to zero with the same steps. At each step, we wait 
long enough to allow the system to reach a stationary state. Then, for the discrete times 
ti = iAt, we measure system characteristics such as the average velocity < v >. 

The symmetries of the equation of motion lead to the existence of nonstationary solutions 
called uniform shding states. They are characterized by the fact that every odd-numbered (or 
even- numbered) atom performs the same motion, but is shifted in time. That is, X2n+i (t) = 
X2n-i{t + Ti) (or X2{n+i){t) = X2n{t + T2)) for n = 1, 2, • • •. It is clear that Ti = T2, 
vi — V2 — V, where vi and V2 are the average velocities of the odd- numbered and the even- 
numbered atoms, respectively. Thus, we need only two dynamical hull functions, fi and /2, 
to describe the motions of all atoms:. 




(2) 



transition by breaking of analyticity occurs [9], Ylf=i sm{2nxj) = 0, so in the case of no 



X2n-l{t) = 'ipo+'ip + (2n - l)uj + vt + 



fi[iP + {2n - l)u; + vt], 



(3) 
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X2n{t) ^ i> + ^nU + Vt + 



f2{ip + 2ruj + vt), 



(4) 



where ipo + ip + (2n — l)uj + vt and ijj + 2nu) + vt are the positions of the {2n — l)th and 
(2n)th atoms for the case of K = 0. The dynamical hull functions /i and /2 describe the 
modulations of the external potential to the atomic positions, ip is an arbitrary phase, and 
ipo is a constant phase. Because of the discrete translation symmetry of Eq. (1), the hull 
functions are periodic, and their periods may be any positive integer rii, i.e. 



It should be noted that fi{ip2n-i + ni -I) ^ fi{V2n-i) and f2{V2n + ni -I) ^ f2{<P2n) when 
we say the period is rii. That is to say, uniform sliding states have infinite kinds of modes. 
This fact is of key importance in obtaining the complete resonance condition. 

In the frame of the mass center which moves with the average sliding velocity v, the 
external potential leads to a time-periodic force acting on the atoms. The frequency of 
this force, the so-called "washboard frequency," is given by the velocity of the mass center 
divided by the period of the potential; i.e., = 2Trv . A resonance may occur if the 
washboard frequency and the eigenfrequencies D,p± of the phonon satisfy the condition 
kiQp+ = k2^w or kiQp- = k2fiw, where ki and k2 are positive coprime integers. To see 
this and to find the eigenfrequencis flp± of the phonon, we (i) use the linear approximations 
sin{27r[(/72n-i + fi{^2n-i)]} ~ sin(27r</72n-i) + 27rcos(27r</72n-i)/i(</?2n-i) and Sin{27r[(/72n + 
/2(v^2n)]} ~ sin(27rv92n) + 27r cos(27r(^2n)/2(v^2n)) (ii) introduce the mean-field treatment, i.e., 
the term cos(27r(^j) is replaced by an averaged quantity 



(ni) consider a small vibration ej{t) of the atoms around their positions, and (iv) use the 
Fourier expansion 



h{.^2n-l + ni) = fl{(p2n-l), 



(5) 



f2{V2n + ni) = /2(</?2n)- 



(6) 




(7) 
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/l(^) 



^ C„2 exp(i 



27rn2 



(8) 



/2M 



(9) 



Finally, we obtain the resonance velocity 



v± = 




{ — [(2 + i^/?)(l + m)±7]} 



(10) 



where {2|;^[(2 + + m) ± 7]}V2 = Q^^, 7 = [(2 + X/?)2(m - 1)^ + 16mcos2(^a;)]V2, 



and 2t:v± = fl^j- It is clear that ki/k2 = ni/n2- This is a complete resonance velocity 
spectrum which relates only to the winding number, i.e., the commensurability, implies that 
the resonance is not a finite-size effect, but a discreteness effect. Prom Eq. (10), it is clear 
that the resonance behavior of the present model is very different from that of the standard 
driven FK model. For a given integer pair {ki,k2), the present model has two possible 
resonance states. For a given resonance state, the value of the average velocity v is fixed, 
so if the value of ki is fixed, then k2 should have two possible values which correspond to 
the two cases where the "+" and "— " signs are taken. So one resonance state should be 
characterized by two integer pairs, {ki, k2) and {ki, k'.^. That is to say, there are two kinds 
of resonance modes in one state. 

When m = 1, we can reduce Eq. (10) to a simpler form: 



where v.. is the real resonance velocity while Vj^ is a fake one. From Eqs. (10) and (11), 
we can clearly find that the external potential not only leads to a periodic driving force but 
also leads to a modulation to the eigenfrequencies. In the high-velocity region, the standard 
FK model shows a behavior similar to that of the diatomic FK model. It is clear that 




(11) 




/3/(47r2) = lim: 




Jo E'{p)dt is the average external potential energy experienced 
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by an atom. In the following, we will see that the dependence of E'^ on F helps us to 
understand better the resonance steps . 

Figure 1 shows the response of the average velocity v of the system to the driving force 
F, where u; = 3/8, 77 = 0.1, X = 1, < F < 0.1, AF = 0.001, and the integral step 
At — 0.01. The hne with sohd squares and the dotted squares without a hne are for the 
case of m = 2. The line with solid circles and the dotted circles without a line are for the 
standard FK model which is plotted for comparision. It is clear that for the diatomic FK 
model, each resonance state is characterized by two pairs of integers, {ki,k2) and {ki,k2), 
which are for small and large atoms, respectively. The main property of a resonance step 
is characterized by one kind of resonance state, but some other resonance states may occur 
at the edges of the step. The ratio ki/k2 is always smaller than ki/k'2, which implies that 
the frequency of the small atom is larger than that of the large one. Figure. 2 shows an 
example of the corresponding dynamical behavior in a resonance state for the case of m = 2, 
where F — 0.015. The thin line is for small stoms and the thick line is for large ones. The 
high peaks are due to hopping from one well to another well, and the low peaks are due to 
neighboring-particle hopping. In this case, the time between two neighboring high peaks for 
the thin line (or the thick line) equals the period of the washboard, and there are four small 
oscillations for the thin line and three small oscillations for the thick line, so the resonance 
state is characterized by (1,4), (1,3). 

Numerical results show that in the high- velocity region, the linear response of ^; to F is 
often punctuated by nonlinear response regions. Figure. 3 shows an example for the response 
of V to F in the high velocity region, where only the process of increasing F is shown: 171 — 2, 
CO = 1/8, rj = 0.1, X = 1, AF = 0.002, and the integral step At = 0.1. It is clear that the 
linear response of f to F is punctuated by resonance steps which result from sub-harmonic 
resonances ( ki > k2). The distances between two successive nonlinear response regions are 
nearly equal. A nonlinear response region is usually composed of one large step and several 
smaller ones. The main resonance properties are characterized by the large step. There 
are two kinds of nonlinear response regions in the high-velocity region. The two kinds of 
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nonlinear response regions occur alternatively. 

To understand better the resonant steps in the high-velocity region, in Fig. 4 we show the 
average external potential energy E'^ experienced by an atom. The values of the parameters 
are same as those used for Fig. 3. From Figs. 4(a) and 4(b), we can easily find that E'^ 
usually oscillates and has different properties for the two kinds of nonlinear responses, which 
correspond to two kinds of dynamical configurations of the moving chain. In one kind of 
the nonlinear response, the time-averaged value of E'^ is higher than that in its vicinity. In 
the other kind of nonlinear response, the time-averaged value of E'^ is lower than that in 
its vicinity. Other interesting phenomena are the facts that the average kinetic energy Ek, 
the average interparticle potential Ep, and E'^ are all constants in the cases of F = 1.5, 3.0, 
4.5, 6.0, • • •, which implies that the moving chain does not feel the existence of the external 
potential in these cases. The results given in this letter also can be used to describe an 
dc-biased one-dimensional Josephson-junction array. 
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Fig. 1. Response of the average velocity v of the system to the driving force F for 
increasing F, where u = 3/8, rj = 0.1, X = 1, < F < 0.1, AF = 0.001, and the integral 
step At = 0.01. 

Fig. 2. Dynamical behavior in a resonance state, where m — 2, F — 0.015. The thin 
line is for small stoms, and the thick line is for large ones. The high peaks are due to hopping 
from one well to another well, and the low peaks are due to neighboring-particle hopping. In 
this case, the time between two neighboring high peaks for the thin line (or the thick line) 
equals the period of the washboard. There are four small oscillations for the thin line and 
three small oscillations for the thick line, so the resonance state is characterized by (1,4), 
(1,3). 

Fig. 3. Response of the average velocity v of the system to the driving force in the high 
velocity region, where only the process of increasing F is shown, m = 2, uj = 1/8, rj = 0.1, 
X = 1, AF = 0.002, and the integral step At = 0.1. 

Fig. 4. Average external potential energy E'^ experienced by an atom. The values 
denoted in the figure are for the external force. The values of the parameters aresame as 
those used for Fig. 3. 
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